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THE GENERALIZED EULER-MASCHERONI CONSTANTS 
INTRODUCTION 


The generalized Euler-Mascheroni constants are defined by 


lim V'' ln k ln +1 M _ _ n i o 

T n m-*» 2^ k - n+ i ’ n °’ 1,2, " 

k=l 


and are the coefficients of the Laurent expansion of 


1 °° (-l) n 7 n (z-D n 

f(z) = -r + E 

2-1 n=0 


They were first defined by Stieltjes in 1885, discussed by Stieltjes and Hermite [ 1 ] , and have been 
periodically reinvented over the years [2,3,4] . 

The rate of convergence is painfully slow (Table 1) so one is forced to seek some method to speed up 
convergence, and of course the most common way is by employing the Euler-Maclaurin formula. Following 
a suggestion made by Edwards [5] in computing the Riemann Zeta function, the sum portion of 7 n is split 

into two sums so that the first (p — 1 ) terms are directly summed and the p to M sum is approximated by the 
Euler-Maclaurin formula as M tends to infinity. Hence, 


P-1 

Tn = S 

k=l 


l nn k + lim 

lr M - * 00 


rM 

l 


ln n x 


dx+ 1/2 


"ln n M ln n p~ 

^ ®2k 

M + p 

A (2k)! 

k=l 


/ ln n x \ (2k-l) 
\ x / x=M 


/ ln n x\ ( 2k 11 
\ x / x=p 



( 1 ) 


where 0 < 0 < 1 and are the BemoulM numbers. Letting M -* it follows that 



Tn 


1=1 ln n k ln n+1 p 1 ln n p B 2k /ln n+1 x \(2k-l) ^ B 2j ^ / ln n x V 2 J 

S gs k ' p +2 p ~h (2k)! V x /x= P ; x=] 


(2j) 

p+k+0 


n (j) 

The chief difficulty in using this formula is determining (ln n x/x) in a usable form. 
Lemma. 


d> / ln n x \ ln n-J x 


dxJ \ x / xJ +1 k=1 


(y-k In x) 


( 2 ) 


where, after expanding the product, put y k = n!/(n-k)!. 
Proof: Let 

m 

Q m = n (y-k In x) 
k=l 

and set 


d m 

dx m 


ln n x 


ln n ~ m x 

„m+l 


<m 


( 3 ) 


and assume that, for some value of m, 


D 1 


,m+l 


ln n x 


= D 


In 


n-m v 


„m+l 


Qm + 


ln n_m x ln m ~ n - l x 


Q' = 

„m+l m 


,m+2 


[n-m - (m+1 ) In x] Q 


m 


+ x In x Qj^] 


and define 


H m = [n-m - (m+1) In x] Q m + X In x Q' 


m 


£ m 


( 4 ) 


Now 
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Qm = 


n i p^nllnx P 2 n!ln 2 x 


(n-m)! (n-m+2)! (n-m+1)! 


+ . . . + (-l) m p m ln m x 


_ y m-l pj Inx + p 2 y m-2 ln 2 x + . . . + (-l) m p m ln IH x 


m „ i„nu 


= y m -y 


( 5 ) 


Since 


(n-m) 


n! 


n! 


(n-m)! (n-m-1)! 


= y m+l 


we obtain, upon substitution of (5) into (4): 


H m = y m+1 * y m [Pi+m+l ) Inx + y m_1 [p 2 + (m+1) p 2 J ln 2 x + ... + (-l) m (-1) (m+1) ln m+1 x 


This is obviously 


m+1 


[y - (m+1 )] Q m = n ( y _k 111 x > = Qm+1 • 
k=l 

Inserting this into (3) proves the theorem by induction. The last term in equation (2) is 


B 2j P^! / b n x \ (2j 

«>' h\ x /x= 


(2j) 

p+k+0 


( 6 ) 


and must now be given bounds. By the Cauchy integral formula, 


ln n (p+k) 

p+k 


(2j) < ( 2!2l 2?rR ln n (p+k+Re 0i ) < (2j)! 2 2j ln n (k+p-2) 


lit R 2j+1 p +k+Re 0i 


P 2j P 


where R = p/2, 


ln n (k+p+p/2 e i0 ) | < 21n n (k+2p) 


and 
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P p+k+p/2 


Inserting this bound into (6) yields 


B 2j (2j)!2 2 J 

(2j)!p 2 i +1 



ln n (2p+k) 

P 




— — — ■ pln n (3p+l ) - : 22Jl ” n ( 3 P' | - 1 > 2 «l)l s 2(2j)!ln”(3p+l) 
p 2 J +2 p 2 J + 1 (2 tt ) 2 J ( 1 -2 1 _2 ->) (pi r) 2 Hl-2 l ~ 2 b 


(7) 


COMPUTATION 


The fourth term of equation (2) requires the construction of three matrices A, B, C and two vectors 
u, v. A is a (2j— 3)X(2j— 2) matrix whose first (L + 1) elements in the Lth row are given by the coefficients of 


L 

H (y-k ln p) 

k=l 


and the remainder set equal to zero. B is a matrix of the same dimension whose Lth row is given by 
n! n! 

( (n-L)! ’ (n-L-1)! ’ ’ ' ' ) 

The elements are, of course, zero if L > n. C is a matrix defined by 


C = W 


for i = 1, 3, 5, 7 ..., and j = 1, 2, 3, 4... 


Next define 


u T = (1 , ln p, ln 2 p, ... , ln 2 J~ 2 p) 


and 




1 

' ’ (2j-2)! 
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With these equations, (2) may be written as: 


> ln n k 
k=l 


ln n+ *p 1 ln n p t* 

— + — — — - v A Cu 

p 2 n — — 


with an error bounded by (7). An APL computer program was written to evaluate 7 n (see Appendix). 

It was observed that the minimum error in computing 7 n occurred in a neighborhood of p=j and a 
sensitivity study indicated that p=j=10 was an optimum choice. Table 2 contains the first 32 Euler- 
Mascheroni constants. 

Table 3 shows that the Laurent expansion provides a very effective means of computing the Riemann 
£ function in a neighborhood of z = 1 . However, the expansion is not a useful method for extending the 
list of zeros of £(z) known today. 

The behavior of the Euler-Mascheroni constants themselves have been the subject of investigation. 
Briggs [6] showed that infinitely many 7 n are negative and infinitely many are positive and Mitrovic 

extended this result by showing that each of these inequalities 7 n < 0, ~t2n-\ < 0. 7 n > *-*> 72n-l ® 

for infinitely many n [7 ] . 

Good [8] recently conjectured that the lengths of the runs of the same sign of A 7 n never decrease. 
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TABLE 1. 7 n VERSUS NUMBER OF TERMS COMPUTED (M) 


M= 10 3 

0.57771558156810 

-0.06936246015836 

0.01416535317172 

7.86463444728723 


M = 10 4 

0.57726566406712 

-0.07235533352434 

-0.00544890007751 

3.31473675482265 


M = 10 5 

0.57722066488224 

-0.07275828112056 

-0.00902762943042 

1.01213466958143 


M = 10 6 

0.57721616479093 

-0.07280893950511 

-0.00959495724544 

0.25241822958923 



TABLE 2. 


n 

0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 


GENERALIZED EULER-MASCHERONI CONSTANTS 


T n 


0.57721566490153 

-0.07281584548368 

-0.00969036319287 

0.00205383442030 

0.00232537006546 

0.00079332381728 

-0.0002387693455 

-0.0005272895671 

-0.0003521233539 

-0.0000343947747 

0.000205332814 

0.000270184439 

0.00016727291 

-0.00002746381 

-0.00020920927 

-0.00028346867 

-0.0001996969 

0.0000262769 

0.0003073682 

0.000503605 

0.000466342 

0.0001044 

-0.0005416 

-0.0012439 

-0.0015885 

-0.0010746 

0.000657 

0.003477 

0.006399 

0.00737 

0.00355 

-0.00752 


TABLE 3. 


X 

iy 

Re£(z) 

Imf(z) 

0.5 

0 

-1.46035450881 

0.00000000000 

0.01 

0 

-0.50929071404 

0.00000000000 

2 

0 

1.64493406685 

0.00000000000 

3 

0 

1.20205690316 

0.00000000000 

4 

0 

1.08232323371 

0.00000000000 

5 

0 

1.03692775514 

0.00000000000 

6 

0 

1.01734306198 

0.00000000000 

0.5 

1 

0.14393642708 

-0.72209974353 

0.5 

3 

0.53273667097 

-0.07889651343 

0.5 

5 

0.70181237117 

0.23103800839 
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APPENDIX 


VEULERlQ ] 7 

7 R+NK EULER MiLiLMiUTiKiN 

[1] LAML«-®.V)* 1 + i2xX*M[2];.WJ:[1] 

[2] 7«-(+/( \M)-{ (L*N + 1 )*A?+1 ) + 0. 5x(L*N)iM 

[ 3 ] T- t -(iM*l + \~l+2*K)x(.L*N-\~l + 2xK)x((,GERM 2 *K) * ( 2*K) GENMN N) + .*LM 

[4] R+N,Y- + \ ( 2+BERNOULLI A+2 ) * r[ A + ~1 +2 x i A 3 * 2 * i A 

[5] a .VA[l3 IS THE N'TH EULER CONSTANT 

[6] a VA [ 2 3 IS THE NUMBER OF TERMS IN THE BERNOULLI SUM 

[7] p M IS THE NUMBER OF TERMS IN THE FINITE SUM 
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NGENMl ij] 7 
V M+-GENM KiVihT 

[13 AM (A-l ) ,K) pO ;I«-2 

[2] 7«-(2,A)p(Apl) ,-iA 

[3] A#[l;3«-n ;l3,(A-2)pO 

[4] LiT+Vl ;I3» . x(«[ j-i ; 3*0 )/tf[I-l ; ] 

[5] 2M?[1;3 ,0) ,[0. 5]0,2 , [2;3 

[6] ¥[I;]*(+/D,(M+l)pO 

[7] +U>I+-I+1)/L 

[8] a GENERATES THE COEF. OF THE DERIVATIVE OF ( LOG X)*N/X 
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7GENMN CL33 7 
7 X+K GENMN N T 

[1] AM (A-l ) ,A) oO ;J<-2 

[2] AfCl;3«-U,l),(K-2)pO 

[3] L:MiI; 3«-(I A'/? TV) ,~1+A/[X-1;3 

[4] -*(A> IM+1 )/L 

[5] a GENERATES MATRIX OF N* (N- 1 ) * . . . * ( N- K+l ) 

[63 A VITA FUNCTION EULER 

7 



^BERNOULLI [Q]7 
7 B+-BERNOU LLI S 

[I] +{N>22)/L 

[ 2 ] 8«-2 2 p 0 

[3] B[ 1 

[4] B[ 2 ]-*-~l *2 

[5] B[ 3 ]-*-l *6 

[6] B[ 4 ]«-~l *3 0 

[7] B[5>1*42 

[8] B[ 6 ]«- - l *30 

[9] B[ 7 ]«-5 *66 

[10] B[ 8 ]■«-"* 6 91*2730 

[II] B[9]«-7*6 

[12] B[ 10]«-~3617*510 

[13] B[ 1 1 ]-*-4 3 867*798 

[14] B[12> 174611*330 

[15] B[ 13]-«-854513*138 

[16] B[ 14 ]«- 236364091*2730 

[17] B[ 15]«-8553103*6 

[18] B[ 16]«- 2 . 374946510*870 

[19] B[ 17 ]«-8 . 615841512*14322 

[20] B[ 1 8 ]«- 7.709321512*510 

[21] B[ 1 9 ]«-2 . 577688512*6 

[22] B[ 2 0 ]«- 2.631527519*1919190 

[23] B[21]«-2. 929994515*6 

[24] B[ 2 2 ]«- - l . 929658516 

[25] ->0 ;B«-B*B 

[26] L: 'N MUST BE LESS THAN OR EQUAL TO 22’ 

[27] n FIRST 22 BERNOULLI NUMBERS 
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